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I. INTRODUCTION 

The microscopic foundations of thermodynamics are well established using the Gibbs hypothesis of statistical en- 
sembles maximizing the Shannon entropy [1] .When the thermodynamic limit can be taken, the various Gibbs ensembles 
for infinite systems converge to a unique thermodynamic equilibrium [28]. However, many systems studied in physics 
do not correspond to this mathematical limit of infinite systems [13] and, in fact, finite systems are now, per sc, a 
subject of a very intense research activity, from metallic clustcrs[3, 4] to Bose condensates [5, 6], from nanoscopic 
systems[7] to atomic nuclei [8, 9] and elementary particles[10]. The question thus arises: can the equilibrium of a 
finite systems be defined. 

A priori, the Gibbs concept of statistical ensembles of replicas which is applicable for an arbitrary number of 
particles, seems an ideal tool to define the thermodynamics of finite systems. However, in a finite system the various 
Gibbs ensembles are not equivalent [11] and lead to different equilibria which physical meaning and relevance has to 
be investigated. In the different physical cases, which in the following encompass both the case of isolated systems 
and of systems in contact with a finite or infinite reservoir [29], the identification of the relevant statistical ensemble 
is a key issue. 

From a macroscopic point of view, a common interpretation of a statistical ensemble is an infinite collection of 
infinite subsystems of the studied infinite system in a specific thermodynamic situation. The required independence 
of the different subsystems is insured by the thermodynamic limit. Within this interpretation, a single system can 
be considered as a statistical ensemble and thus can be disciissed in terms of equilibrium. This kind of equilibrium 
is not relevant for a finite system since i) the interface interactions between subsystem cannot be neglected, ii) the 
procedure of coarse-graining modifies the entropic properties of the system and iii) a finite system does not lead to 
an infinite ensemble of subsystems. In fact a single realization of a finite system cannot be discussed in statistical 
physics terms. 

An alternative viewpoint is given by the Boltzmann ergodic assumption. In this interpretation the statistical 
ensemble represents the collection of successive snapshots of a physical system evolving in time. The equivalence 
between this time average and the Gibbs ensemble is then insured by the ergodic theorem [30]. This interpretation 
however suffers from important drawbacks. First, not only a proof of the ergodic hypothesis under fairly general 
conditions is lacking, but even for a truly ergodic Hamiltonian, a finite time experiment may very well achieve 
ergodicity only on a subspace of the total accessible phase space [14]. Moreover, ergodicity applies to confined systems 
and thus it requires the definition of boundary conditions when the thermodynamic limit docs not apply. Then the 
statistical ensemble in general explicitly depends on the boundary conditions and we will discuss in this article that 
an exact knowledge of the boundary corresponds to an infinite information and is therefore hardly compatible with 
the very principles of statistical mechanics, as a reduction of the many body information to a (small) number of state 
variables. 

Finally, in many physical situation the ergodicity ideas do not apply. Indeed, many physics experiments do not 
follow the time evolution of a single system, but rather concern averages over a great number of events, i.e. of physical 
replicas of systems experimentally prepared or sorted in similar way which are then observed at a given time. In such 
a context there is a priori no connection between the measuring time and the time it takes for an ergodic system to 
visit evenly the energy shell. Moreover, the systems experimentally accessible are often not confined but freely evolve 
in the vacuum, as it is notably the case for atomic and nuclear clusters and high energy heavy ion collisions. These 
open transient finite systems are isolated but never stationary. The concept of a well defined stationary equilibrium, 
uniquely defined by the variables conserved by the dynamics in an hypothetical constraining box, is certainly not 
useful for these finite systems. 

However statistical approaches, expressing the reduction of the available information to a limited number of collective 
observables, can still be of some pertinence for such complex many body systems [15]. The physical meaning of such 
statistical ensembles is different from the case of an ergodic system or a macroscopic coarse graining. The standard 
assumption is that the dynamics is sufficiently complex (chaotic or mixing) [16-18] such that, repeating an experiment 
many times, the ensemble of events dynamically populates a phase space in a "democratic" enough way such that the 
gross fcatiires of the ensemble of events arc dominated by few collective variables. In this case, the maximum entropy 
postulate cannot be justified from the ergodic theorem but has to be interpreted as a minimum information postulate 
which finds its justification in the complexity of the dynamics independent of any time scale[l, 15]. 

This informaticm theory approach is a very powerful extension of the classical Gibbs equilibrium: any arbitrary 
observable including time odd quantities can act as a state variable, and all statistical quantities as equations of state 
and phase diagrams can be unambiguously defined for any number of particles [19]. The price to be paid for such a 
generalization is that the density matrix continuously evolve in time as soon as the constraining observables are not 
restricted to conserved quantities, meaning that we have to take into account the time dependence of the process. 
A well known example of application of the information theory to a dynamical process is given by time dependent 
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mean field theories, whieh can be viewed as the solution of a maximum entropy variational problem under the time 
dependent constraint of the knowledge of all one body observables [15]. 

In this article, we develop an information theory approach to the physics of finite systems evolving in time. We 
show that such an approach provides a description of the thermalization process. In presence of non linear dynamics 
it may lead to deviation from the Boltzmann-like exponential distributions and thus might provide an explanation 
for the appearance of non-extensive statistics [2]. In a mean-field context this dynamical extension envolves the RPA 
matrix leading to the introduction the collective degrees of freedom. For unconfined systems, we will show that their 
finite size at a finite time can be accounted for by introducing, in the Shannon information, additional constraints 
describing the system's compactness (size and shape). This provides a natural solution to the boundary condition 
problem. The same formalism describes incompletely known boundary conditions, confinement by a potential, and 
unconfined self bound or unbound systems. The time dependence of the process naturally leads to the appearance of 
new time-odd constraints or collective flows. In the case of an ideal gas of particles or clusters we will show that the 
extension and flow constraints forms a closed algebra together with the energy operator allowing an exact description 
of the free expansion of the system in vacuum. 

II. STATISTICAL EQUILIBRIA 

Let us first recall the standard formalism of statistical equilibria which we will then extend to the case of time 
dependent processes involving finite systems. 

In order to describe a statistical ensemble, we introduce either the classical density D {q, in the many-body 

phase space associated with all the particle positions Q = {qi} and momenta P = {pi} or the quantum density matrix 
= X)(n) 1*^"^) P^"'^ <;|\l/(") I where |\1>(")) are the state of the different events (n) and p^"^ the associated occurrence 
probability. According to the Gibbs hypothesis, equilibria are maxima of the Shannon entropy [1] 

S= -Tr£»log£», (1) 

where Tr means an integral over the coarse-grained classical phase space, (^Q, , or the trace over the quantum Fock 
or Hilbert space of states 1^1/). In this article we implicitly use units such that the Boltzmann constant k = 1. 

When the system is characterized by L observables, A = {A^}, known in average < A^ >= TiDA^, one should 
maximize the constrained entropy 

S' = S-J2^i<^i> 

e 

where the A = {A^} are L Lagrange multipliers associated with the L constraints < Ag >. The Gibbs equilibrium is 
then given by 

D^ = ^exp-\1, (2) 

where X.A = ^i^e and where is the associated partition sum insuring the normalization of D^. 

To interpret the Gibbs ensemble as resulting from the contact with a reservoir or to guarantee the stationaritj of 
the equilibrium eq.(2), it is often assumed that the observables Ag are conserved quantities such as the energy H or 
the particle (or charge) numbers iVj. However, there is no formal reason to limit the state variables to constants of 
the motion and, in fact, the introduction of non conserved quantities is a way to take into account some non ergodic 
aspects. Indeed, an additional constraint reduces the entropy, limiting the phase space and modifying the event 
distribution. This point will be developed at length in the next sections. 

It should be noticed that the formalism recalled above encompass microcanonical thermodynamics[8] which can be 
obtained from the variation of the Shannon entropy eq.(l) in a fixed energy subspace with no external constraints. 
In this case the maximum of the Shannon entropy can be identified with the Boltzmann entropy 

max{S)=\ogW{E), 

where W is the total state density with the energy E. 
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The microcanonical case can also be seen as a particular Gibbs equilibrium (2) for which both the energy and 
its fluctuation are constrained. This so called Gaussian ensemble in fact interpolates between the canonical and 
microcanonical ensemble depending upon the constraint on the energy fluctuation [12]. The same procedure can be 
applied to any conservation law so that the Gibbs formulation (2) can be considered as the most general statistical 
ensemble. 



III. MULTIPLE TIME STATISTICAL ENSEMBLES 



Eq. (2) represents the standard statistical description of a physical system at a given time. Indeed, if some of 
the observables are not constants of motion, then the statistical ensemble (2) is not stationary, but will evolve in 
time. This specific role played by time stresses that this information theory "equilibrium" cannot be justified using 
the usual ergodic arguments, but through the fact that the global features of the replicas of the considered system 
are characterized, at a given time, by few observables, i.e. that the information is concentrated in few degrees of 
freedom. In many physical cases one can clearly identify a specific time (" freeze out" time) at which the information 
concentrated in a given observable is frozen (i.e. the observable expectation value ceases to evolve or presents a trivial 
dynamics). However this freeze out time may be fluctuating or different for different observables. For example for the 
ultra-relativistic heavy ion reactions two freeze-out times are discussed, one for the chemistry and one for the thermal 
agitation. To solve these questions we need to introduce time as an explicit variable and define a statistical ensemble 
constrained by informations coming from different times. 



A. Formulation of the issue 



Let us assume that the evolution of an ensemble can be written as 

dtD = J'[D], (3) 

where T[D] is a functional of the density matrix D. This is a very general dynamical evolution since it includes 
classical and quantal Hamiltonian evolutions 

dtD^{H,D}, (4) 

where H is the system Hamiltonian and {., .} are Poisson bracket in classical physics and commutators divided by 
ih in quantum physics. Eq. (3) also includes non-linear approaches such as mean-field approximations and more 
generally variational treatments for which H is replaced by an effective operator which depends upon the actual state 
D : H ^ H[D]. Eq. (3) also includes the stochastic extensions of such approaches, D being the ensemble average of 
the stochastic evolutions. 

Let us now suppose that the different informations on the system, < Ag >, are known at different times, tf. 

<Ae>t,= TrD (i^)i^ 

A generalization of the Gibbs idea would be that at a time t the least biased state of the system corresponds to the 
maximum of the Shannon entropy, considering all informations as constraints. Causality arguments imply that this 
time t should be larger or equal to all the tg. It should be noticed that in the case of a Hamiltonian evolution (4) 
because the entropy is a constant of motion this remark has no implications and the maximization can be performed 
at any time leading to the very same result. In the case of dissipative systems the entropy grows and eventually 
saturates. In this case the Maximum Entropy principle has to be applied at the time of entropy saturation (or the 
observation time if it occurs before) . 

The maximization of the entropy at time t with the various constraints < A^ known at former times t(, 
corresponds to the free maximization of 

L 

S' = S{t)-Y,^t<A,>t, 

i=\ 

S' = -TV (^D (i) log D{t) + J2 ^'^'^ (^^)) ' (5) 
where the Xi are the Lagrange parameters associated with all the constraints. 
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B. Minimum information under fluctuating-time constraints 



Let us first assume that the various times only slightly differ by Sti = t — ti, thus we can use the equation of motion 
(3) to link the various times by 

D (t,) = D{t)~ St,T[D (t)] 
in order to explicitly write the constrained entropy as a function of a unique density D = D {t) 

S' = -Tr ^log ^ + ^^^^ - 5U\iAeT[b\j . (6) 
Computing the variation 5S' of S' induced by a modification of the density matrix D ^ D + 6D we get 

5S' = -TrSD I log£) + 1 + ^ XgAe + ^ Xi5teBe[D]A , 

where the operator are related to Ai and to the functional derivative A^[i^2][-D] = 9^[i][-D]/9£)pj by 

B[i][D] = -Tr[2]i[2]A^[2,i][^], (7) 

where the indices [i] indicate the space over which the operators are acting and where Z)^j is the transposed matrix 
when it applies, i.e. in quantum mechanics. In this case, if we introduce a base {|/)} to explicit the trace, eq. (7) 
reads: 

Bu[D] = -Y.^K.'-^^. (8) 

The minimum biased density matrix (solution of SS' = 0) is given by 

1 ^ 

(*) = ^77Y - E ^^^^[^A (*)]. (9) 
where we have introduced a modified observable 

A',[D^^{t)] = [Ae + 6teMDxm) (10) 

which takes into account the time difference between the various observations. The associated Lagrange parameters 
are defined by the equations of states 

< A, >t,= Tr (0 - SteJ^iDx (*)]) A, (11) 

The partition sum is defined as 

L 

it) = Tr exp - ^ A.i^l^x (12) 

The equation < A'g >t= —dLogZx{t) / dXi always holds but a similar relation between the equation of state (11) and the 
partitionsum (12) is valid only if jr[_C)] is a linear functional of D. Insuchacasc < A^ >ti=< A'g, >t——dLogZy^{t)/d\(,. 

The minimum information density matrix eq.(9) looks like a standard Gibbs equilibrium but this formal analogy 
hides important differences. Indeed the constraining observables are modified according to (10) which contains new 
operators Be,. If the original observables A(, were time even, the operators B^ are time odd, and therefore correspond 
to non stationary situations. 
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This can be more easily seen if we interpret eq. (9) as an extended Gibbs equilibrium under the set of constraints 

^ L L 

where the parameters vi contains the time information V)- = Stf Xf . 

Because of the possible D dependence of the Be operators, the distribution {t) might deviate significantly from 
the usual exponential behavior. This opens an interesting possibility to encounter non Gibbsian statistics. In the 
literature [2], non Gibbsian information kernels arc usually derived from a modification of the entropy as, for example, 
in the case of a Tsallis distribution. Here these anomalous statistics might be obtained from the time dependence of 
the studied system. In the following sections we will elaborate more on this subject. 



C. Hamiltonian evolution 



Let us illustrate the above results in the case of a Hamiltonian evolution (4). Then the entropy is a constant of the 
motion S = -Tr (^{H, D {t)}{logD + 1)) = Tr (^H{{logD + 1), D {t)}^ = so that the minimum biased trajectory is 
independent of the time t at which the entropy is maximized. Introducing 

and using the cyclic invariance of the trace (or the by part integration in the phase space integral for the classical 
case), eq.(6) can be written as 



S' = -Tr (^D (^ogD + ^ Xde + ^ 6UXi{H, j j 



(14) 



In Eq. (14) we can introduce 

in agreement with the Heisenberg picture, in which the observables would evolve from time te up to a common time 
t. The generalized information theory result, i.e. the extremum of the variation 5", is given by 

^A(*) = ^exp-A.i'. (15) 

In this case the A' operators do not depend on D so that the distribution (15) is similar to a standard Gibbs 
equilibrium. The Lagrange parameter are defined by the constraints 

< Ae >t,= TrD {U) A, = TtD-^ {t) A', = -dx, log (t) . 

Eq. (15) can also be interpreted as the introduction of additional constraints 

Be = {H,Ae} 

and additional Lagrange parameter associated with the out of equilibrium minimum biased density matrix 

^A> = ^exp-A.i-i7.1 (16) 

The equations of state are given by < Ag >t= —dxg logZ^ ^ {t) and < B^ >^= —d^^ \ogZ-^ ^ {t) . Then the specific 
case (15) is obtained requiring V(, = XeSte, the A's being defined by the constraints < Ae >j^=< Ae >t +Ste < Be >f 
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D. Minimum information under multiple-time constraints 



This procedure can be easily extended to longer time intervals. To do that, we have to integrate the evolution from 
te to t 



D{te)=D{t)- [ dt'T[D{t')]. 
The variations of the density matrix at various times are related by 

si) {t') = 6D (t) - ^ dt"j^[D (r) + si) (r)] - j'liD (r)], 



It- 

which leads to the relation 



SD[^{t')=SD[^it)- rdrTr[2i^5H^i!!)l5i) (D, (17) 
Jt' dDI^^ 



where the indices [i] indicates the space over which the operators are acting (see equation (7) ). This equation can 
be integrated giving 



<5^[i] it') = <5^[i] it) - J2 ^— r^Trp] A^fij {t, t') cJ^p] (i) , 
p=i 

where the matrices Aljf 2] defined iteratively by 
with Al[°^2] = 1- If 

we introduce this expression into the constrained entropy S' eq.(5) to compute the variation 5S' 

,we get 

5S' = -TtSi) it) (^ogDit) + 1 + E + E ^^—^B^i'^ ^^)) ) ' 
where the operator B^^^ are related to and to the matrices jj'j] through 

(i> ti) = -T^n2]Ae^,M[2,i] • (18) 

In the special case of a Hamiltonian evolution, the operators B^f^ are simply p— uple commutators M^^^ = {H, B^p~^^} 
with B^"^ = A. Therefore in this special case they do not depend neither on time nor on D. Coming back to the 
general case, as in eq.(lO) above, we can introduce the new observable 



A',{t,te) = A, + f2^-L-^B^^^ {t,U) 
The minimum biased density matrix is given by 



p=l 



1 



(*) = exp -X.A' {t, U) . (19) 



where the partition sum [t) = Trexp— A.A' {t,t(i) insures the normalization of {t). The Lagrange parameter 
are defined by the constraints 



< Ae >t= Tr (^t)^ (t) - ^* dt' J^[D^ {t')]^ i,. 



(20) 
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which only for a linear dependence of diD on D such as in the particular case of a Hamiltonian evolution give 
back the standard expression for the equations of state < Ai >ti= —d\f.\ogZ^(t) since in this case < >ti = 

Tr£)j^ [t) A'^ {t, te) =< A'^ >t. It is interesting to note that if the evolution is Hamiltonian, the A'^ represent the time 
evolution of the constraining observables A( in the Heisenberg representation A'^ {t, te) = A'^ (Ate) = e~'^^*'^ Age^^*^^ 
where Ati = t — te. 

Coming back to the general case, as for eq.(9) above, eq. (15) can also be interpreted as the introduction of 
additional constraints B^^^ and additional Lagrange parameters v^^^ associated with the time evolution of the system 

1 oo ^ 

%- = ^e^p-^-^-E^^'^-^^''^ (21) 

where the Lagrange parameter v^^'' arc related to Xi by i'^^' = ^^^—^r—^i ■ 

It is important to notice that Eqs. (19) or (21) provide exact solution of the complete many body evolution 
problem eq.(3) with a minimum information hypothesis on the final time t having made a set of observations < Ae > 
at previous times ti. We can see from cq.(21) that in general an infinite amoimt of information, i.e. an infinite 
number of Lagrange multipliers are needed if we want to follow the evolution of the density matrix for a long time, or 
if the time interval between each constraint is long. However, if we are only interested in short time scales the series 
will rapidly converge. Moreover in the next sections we will show that different interesting physical situations exist, 
for which the series can be analytically summed up. In this case, a limited information (the knowledge of a small 
number of average observables) will be sufficient to describe the whole density matrix at any time, under the unique 
hypothesis that the information was finite at a given time. 

In the next sections we will illustrate this theory with two representative examples: a dissipative non-Hamiltonian 
dynamics and a self consistent case. 



E. Application: Brownian motion 



The description of the dynamics of quantum systems interacting with their environment is, in general, a very 
difficult task due to the system-reservoir interaction that usually involves a huge or even an infinite number of degrees 
of freedom. The Liouville-von Neumann equation for the total closed system is therefore useless for a reasonable 
description. The standard prescription is then to construct effective equations of motion for the reduced density 

matrix of the system by tracing over the environmental variables in the exact dynamics. The resulting quantum 
master equation includes dissipative and stochastic terms that take into account the irrelevant degrees of freedom in 
an approximate way. 

One prototype of system-reservoir models is the Caldcira-Lcggett model [20] for a Brownian particle of mass to, 
with coordinate f and momentum p , in a bath consisting of a large number of harmonic oscillators. After performing 
a series of approximations [20], the density matrix of the Brownian particle is found to satisiy the following Liouvillian 
quantum master equation : 



dtt) = J^[D] = 



r ^2 



ih 



2m' 



D 



p,D 



d_ 



r,D 



(22) 



The first term describes the free Hamiltonian dynamics. The second term, proportional to the friction coefficient 

7, is a dissipative term associated to the mean coupling to the bath. The last term, with the diffusion coefficient d 
satisfying the Einstein relation d = 2to7//3 where /?~^ is the bath temperature, describes thermal fluctuations. 

Let us apply the general formalism of section III D to a one dimensional Brownian particle prepared, at a time , 
with a mean kinetic energy < K >to=< §^ >*(,, and which is observed at a later time t. Following the procedure of 
section HID, the information theory ansatz for the density matrix D(t) is a generalized statistical equilibrium with 



the constraining observables 5 
easily show that 



K and ^[fj^ 



-TrroiB, 



(p-i) 



[2] ^[2] 



9^[2]/9D^y With the irreversible dynamics (22), we 



^(1) = 47X - - 

TO 

Bip) = 47B(f-i) = (47)^-' 



(23) 
(24) 
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Using (19), the density matrix D{t) can be written as 

where Z'^ {t) = Trexp— A^' {t) and where the modified constraining operator K' {t) can be actually resumed as 

K'U) = y (^-^o)^ ^(p) = e47(t-*o)^ _ J_ ( e47(t-to) _ . (25) 

The Lagrange parameter A fulfills the equation of state < K' >t= —dx log Z'^ {t) which gives after an explicit calcu- 
lation of the partition sum 

1. = ^4i(t-to) <K>t. (26) 

Since the evolution Kernel T[D] is a linear functional of D, the relation < K' >t=< K >tgholds (see eq.(20)). Taking 
advantage of the explicit expression (25) of the observable K'{t) we get 

< k >t=< k >t, e-^T(*-*°) + ^ (l - e-4T(*-*°)) (27) 

This expression corresponds to the exact mean value of the kinetic term at time t, as deduced from the equation of 
motion of < k > induced by the quantum master equation (22 ) 

d<k> ^ . d 

; = -47 <K > -\ 

at m 

which clearly shows the damping role played by 7 and the fluctuation ensure by d. 

Let us now interpret the above results. Since k' {t) is a combination of the unit and kinetic operator, the density 
matrix can be recasted as 

Dx(t) {t) = exp -X{t)i^ {t) k, 

with //(f) = exp{4"/{t — to)). In the above equation we have make explicit the fact that the Lagrange multiplier A 
is time dependent as can be seen solving Eq. (26) with the help of Eq. 27. This density matrix can be interpreted 
as a canonical equilibrium with a time dependent temperature defined by {t) — {t). Indeed the equation of 
state defining \{t) implies T {t) = 2 < k >t which is the standard result for an ideal ID gas. Using eq.(27) we can 
work out the time dependence of the temperature 



T{t) exponentially relaxes from the initial temperature 2 < k >tgto the bath temperature with the correct 
characteristic time I/47, in agreement with the exact dynamics[20]. 

Thus applying the generalized Gibbs ensembles for time dependent processes to the problem of a quantum particle 
in contact with a reservoir we have been able to describe the general evolution as a succession of canonical statistical 
ensembles correctly relaxing toward the expected equilibrium. This excmplc illustrates the potential of the developed 
formalism to statistically describe time dependent out-of-equilibrium phenomena. 



F. Application: Mean-field and non-linear dynamics 



As an illustration of a non-linear dynamics, let us study the case of self-consistent approaches which are of the mean- 
field type. To be more specific, we consider a quantum system characterized by its one body density p , which can 
either be seen as the density matrix projected over all particles but one (/5[i] = Tr[2^..._Ar] ^[1,2,. ..,Ar]) or the expectation 
value of a generic one body operator pj^j =< ajUj > where 0^(07) are creation (anihilation) operators of a particle 
in the orbital J {I). Since we are studying a one-body approach we will restrict the discussion to one-body operators 
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The time dependent mean-field dynamics is given by 

dtp = {W\p\rp}, (28) 

where W[p\ is the self-consistent mean-field Hamiltonian. Using a variational approach[15] W[f)\ can be related to the 
functional derivative of the energy E ~< H > : W[p][i] = dE/dpJ^ i.e. Vl^[/5]7j = dE/dffjj. If wc consider a small 
deviation from an equihbrium solution, p = po + 6p the dynamics of 6p follows the time dependent linear response 
(RPA) equation[21] 

St%] = W],<5p}[i] + J^l/[i,2][/5](5/5[2],/5l , (29) 



[2] 



[1] 



where V[i,2] = 9W[p]^^/dp'^^, i.e. ViljkIp] = dW[p]ij/dpKL, is interpreted as the residual interaction. It is more 
convenient to introduce the Liouville space, considering the density matrices p as vectors || p » of components 
a = and one body operators A as dual vectors <SC A || using the scalar product ^ A || p ^= TrA+p [31]. In this 
representation the dynamics of 6p can be writtcn[22] 

idt II Sp >= H II 6p >= {W-AfV) II 6p > (30) 

where the RPA matrix H has a commutator structure W jj • >=|| [W, •] >, AT jj ■ »= || [p, ■] > and V jj • »=|| 

It is easy to verify that the matrix A^[i,2] = 9J^[i]/dp^^ introduced in section HID is — i-times the RPA matrix H 

iM.ij,KL = Hij^KL = WikSlj — SilWkj + Vil,mkPmj — PimVml,jk 

while the new oper ators B( introduced to take into account the propagation of the information < Ai> over the time 
interval 5ti = t — ti are given by 

«:BH|=i«^^ II W=i«^H| (W-ATV) (31) 

To get a deeper insight into this relation, let us introduce the energies uj^ and eigenstates 5pv of the RPA matrix [32] 

oJu II 5pu'> = H II Spi, » (32) 
< II W = < C„ II (33) 

where <C || represents the dual basis -C C^ji \\ 5pv 5v'v The operators can be expressed as 

<^Bc\\=i'^uJu^Ac\\5p^-:>-^C^\\ (34) 

Expanding the <C ^4^ || over the same basis 

<^H|=^«;A^ ||(5p, »<CJ| (35) 

we observe that Af, and are akin to conjugated operators. To make this relation more explicit let us assume that the 
operator A(, excites mainly one state vg, with a real frequency with an excitation amplitude <C Ai || 5pv^ a^^- 
Taking care of the fact that RPA eigenstates appears by pairs [21], Ai and B^ can be written as (see appendix C for 

details) 

A(, = a^fi^^+ alfi+^ (36) 
B(, = iiv^^aveC^e - iLOv^alfi^^ (37) 

If Ag is interpreted as a collective coordinate associated with the creation of a collective mode through the (5+ operator, 
then the extra constraint Bi taking care of the time dependence has the structure of the associated momentum. 
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It is interesting to notice that if the RPA matrix does not evolve in time, the above treatment can be iterated to 
propagations over a finite time t. In this case 4C B^,^^^^ \\= i <C B^^^ \\ H so that if we consider a single mode 

Sf^' = {-clfA, (38) 

Sf^+i) = {-^IfB, (39) 

then the time dependence of the process can be accounted through a simple time dependence of the constraining 
operator 

A'^ (t) = Ai coauji,^ (t -ti) + BiLo:^^ amuj^^ (t-ti). 

which represents an oscillation in the RPA collective phase space A^ and B^. 

Thus applying the proposed formalism to non-linear dynamics we have shown that the Gibbs ensembles description 
for time dependent processes can take into account new degrees of freedom. In the mean-field example, the RPA 
modes appear naturally and the generalized ensemble correctly describe their collective oscillation. This demonstrates 
how the developed formalism allows to go beyond the standard approaches to take into account dynamical aspects in 
a statistical manner. 



IV. UNBOUND SYSTEMS AND THE DYNAMICS OF THE EXPANSION 

An interesting application of the developed formalism is given by the case of finite systems in the presence of a 
continuum. For such physical systems the boundary condition problem is far from being trivial as we discuss below. 
We will show that the only consistent way to describe the boundary of unbound systems is to introduce in the Gibbs 
equilibrium at least one observable related to the system's extension (size and shape). Since these observables are 
not constants of the motion, such states are not stationary: in the absence of a confining force, the size of the system 
increases with time and collective flows are generated. An especially relevant physical case is given by nuclear collisions 
[9, 26] which produce a transient correlated high density state that subsequently freely expands in the vacuum. We 
will show that the out of equilibrium expansion can be accounted for at any time within a finite amount of information, 
with the natural emergence of time odd collective flow observables among the state variables. 

A. Boundary condition problem 

The information theory formalism is valid for any system size and thus can be, a priori, applied to finite systems. 
However, as soon as one A^ (e.g. H) contains diff'erential operators such as a kinetic energy, the standard Gibbs equi- 
librium eq. (2) is is a priori not defined, unless boundary conditions are specified. The same is true for micro canonical 
thermodynamics, since this latter can be seen as a special case of the more general Gibbs formalism. Boundary 
conditions arc irrelevant at the thermodynamic limit (when it exists), or if the system is bound: in both cases the 
definition of a boundary corresponds to negligeable surface effects. However, in a finite system in the presence of a 
continuum this is never the case, and in this respect the statistical physics of small unbound systems is ill-defined. 

1. Systems in an external potential 

Physical situations occur in which an external potential is used to confine the system, as for recent studies of bosons 
or fermions in traps. In this case boundary conditions are not a problem since the confinement is directly insured by 
the external potential and the boundaries can be rejected to infinity as is the case of self bound systems. 

However, the external potential directly appears in the system Hamiltonian and thus it is clear that the statistical 
properties of the system directly depend upon the considered potential. For example, systems are often trapped in a 
one body harmonic potential well U = u„ where the sum runs over all the particles n and where u = kr^ jl. The 
corresponding canonical ensemble reads 

t)0M = ^ exp -/3(JJ + fc^V2), (40) 

where W = is the total radius squared. 

We can see that the confinement via a potential is equivalent to an extra constraint defining the compactness of 
the system through the introduction of a Lagrange multiplier given by the oscillator strength. [33] 



12 



2. Systems in a fictitious box 

In the absence of a physical constraining potential, systems presenting states in the continuum are often confined 
through boundary conditions: a fictitious container is introduced in the theoretical description of unbound systems to 
limit the Hilbert space in order to be able to define a statistical ensemble [23] . The introduction of such an unphysical 
box not only directly affects the energy spectrum and thus the thermodynamics properties of the system, but also 
leads to an intrinsic inconsistc!iicy of the statistical theory, as we now show. 

Let us write the eigenvalue problem for an Hamiltonian operator H containing a non-local term such as the kinetic 
energy. To be solved we must specify boundary conditions which directly modify the non-local operators. As an 
example, let us consider the standard case of the annulation of the wavefunction on the surface S of a containing 
box V. Introducing the projector, Pg , over the surface S and its exterior, the boundary conditions reads Ps\^) = 
or, using P| = Pg, {^\ Ps \^) = 0. The eigenvalue equation is then equivalent to a variational principle where the 
boundary and the normalization constraints have to be taken into account through the introduction of two Lagrange 
multipliers bs and E 

6 ((*| H I*) - I*) - 1) - 6s (*| Ps I*)) = 0. (41) 
leading to a modified Schrodinger equation 

H\'^Es)~bsPs\'<i^Es)=E\^Es) (42) 

Equation (42) shows that the energy E and wavefunction \'^es) directly depend upon the boundary conditions. 
Let us now extend this discussion to mixed states. The boundary condition Pg j^f^")) = is exactly equivalent to 

the extra constraint < Ps >= TvDPs = 0. If we note again A the observables characterizing a given equilibrium, the 
density matrix including the boundary condition reads 

^xs = ^ exp - A.i - bsPs (43) 

which shows that the thermodynamics of the system does not only depend on the Lagrange multiplier bs, but on the 
whole surface S. 

For the very same global features such as the same average particle density or energy, we will have as many different 

thermodynamics as boimdary conditions. More important, to specify the density matrix, the projector Ps has to 
be exactly known and this is in fact impossible. The nature of Ps is intrinsically different from the usual global 
observables A^. At variance with the A^, Ps is a many-body operator which does not correspond to any physical 
measurable observable. The knowledge of Pg requires^ the exact knowledge of each point of the boundary surface 
while no or few parameters are sufficient to define the A^. This infinity of points corresponds to an infinite amount of 
information to be known to define the density matrix (43). This requirement is in contradiction with the statistical 
mechanics principle of minimum information. Thus eq.(43) is unphysical. 

This incoherence should be solved by a statistical treatment of our knowledge on the boundary condition. The den- 
sity matrix should be an average over the boundary surfaces of D^^ weighted by the associated boundary probability 
derived applying information theory on our actual knowledge of the boundaries. A simple way to perform this task 
is presented in the next section. 



3. Incomplete knowledge on the boundaries 

One way to get around the difficulties encountered to take into account our incomplete knowledge on the boundaries 
is to introduce a hierarchy of observables describing the size and shape of the matter distribution. 

If the system is unbound, the unconstrained equilibrium corresponds to an infinitely expanded system. However, 
the relevant physical problem in general concerns the state of the system before this asymptotic stage, when it has 
only expanded to a finite volume. This transient stage is constrained by a finite size, i.e. a finite average of an 
observable related to the system's size. 

For example, if only the average system size < P^ > is known, information theory requires the introduction of one 
additional Lagrange multiplier imposing this information, i.e. one of the A^ describing the system is R^. [34] If the 
additional information reduces to the energy, the Lagrange multipliers associated with the state variables < H > and 
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< > arc respectively Xh = (3 = l/T, the inverse of a temperature and A/f2, which has the dimension of a pressure 
when divided by a typical scale Rq and by the temperature, \r2 — (3PRq. The minimum information principle implies 

Dp,p = -J— exp -a (h + PRoR^) , (44) 

which is akin to an isobar canonical ensemble. It is interesting to remark that the term proportional to R^ can be 
equivalently viewed as a mean square radius constraint or as an harmonic external potential modifying the effective 
Hamiltonian of the system as in the equilibrium (40 ). Similarly the expectation value < H + PRoR? > can be 
interpreted as an enthalpy or as a constrained energy. 

A typical application of eq.(44) is given by the freeze-out hypothesis proposed for the unconfined transient finite 
systems produced in atomic or nuclear collisions[9, 26]. At a given time the main evolution (i.e. the main creation 
of entropy) is assumed to stop. Before this freeze-out time the evolution is assumed to be complex enough such that 
all the partitions compatible with common gross features are freely (i.e. with no bias) c;xplored. After the freeze-out 
time, partitions are supposed to be essentially frozen because of the lack of interactions. Typically thermal and 
chemical equilibrium is assumed, meaning that the information on the energetics and particle numbers is limited to 
the observables < H > and < Nj > for the different species / [9, 23]. 

The freeze-out occurs when the system has expanded to a finite size. Then at least one measure of the system's 
compactness should be included among the collective variables characterizing the statistical ensemble. The limited 
knowledge of the system extension leads to a minimum biased density matrix given by eq.(44). [35] 

Such an ensemble does not suffer from the drawbacks of models introducing a constraining box; i.e. both the 
arbitrariness of such a fictitious box and the conceptual problems of the statistical treatment of boundary conditions. 

To conclude this section, we have shown that all the different ways to describe boundary conditions in finite 
unbound systems lead to the introduction of additional constraints in the set of relevant observables describing 
the statistical ensemble. The more information are known on the spatial extension of the system, the more complex 
are the additional observables. In other words, only Gibbs equilibria with compactness constraints are physically 
meaningful when discussing small systems in the presence of continuum states. Since outside the thermodynamic 
limit statistical ensemble are never equivalent, the thermodynamics of finite unbound systems always depends on the 
compactness conditions. 



B. Unconfined finite ideal gas 



Let us first consider the case of a finite ideal gas, i.e. the Hamiltonian is reduced to the one body kinetic term 

^ = — . (45) 



2m 2m 

n 

Let us assume that the first constraining observable is = and so the associated Lagrange multiplier is 
the usual inverse temperature /3. As discussed in subsection IV A, any statistical ensemble is ill-defined unless 

boundary conditions are specified. In the case of an unbound system, one is forced to introduce additional observables 
constraining the spatial extension of the system to a finite size. This is a situation often encountered experimentally: 
a finite system of loosely interacting particles or clusters with a finite extension after an expansion during a limited 
time in an open space. The minimal assumption needed to characterize statistically the ensemble of states, is given 

by one single observable related to the size, e.g. the knowledge at a given time of the mean square radius < R^ > 

(with R^ = J2n^n)- Then we have to introduce the constraining observable A2 = R? associated with a Lagrange 
multiplier Aq in the statistical description. The maximum entropy solution is given by 

6,,, = _i_e.p-,X:(| + ^?„). (46) 

Eq.(46) is akin to a system of non-interacting particles trapped in an harmonic oscillator potential with a string 
constant k = 2Ao//3. From the partition sum, the EOS are easily derived. For example, in the classical case using the 
equipartition and the virial theorem it is easy to derive the system EOS for each particle n = 1, . . . , A'' 

<pi> = ^ ^ ^ a 

2m /3 " 2/3' 
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i.e. 

3rn 3 

Since in the presented derivation, tlic XqR^ term is not an external confining potential but only a finite size constraint 
, the minimum biased distribution (46) is not stationary {{H, D^x^} ^ 0). The system represented by eq.(46) at a 
time to will evolve in time according to the Hamiltonian (45). 

The physical situation we are describing can be explicitly realized experimentally by taking a system in equilibrium 
in an harmonic trapping potential and by suddenly removing the confining potential at a time to. After 

that, the non- interacting particles freely expand ballistically (see section IV A 1). In this case the information theory 
ansatz (46) with Aq = fc/3/2 describes the standard static equilibrium just before the trapping potential is removed, 
and is therefore the exact initial state. This is akin to the actual experimental procedure used in Bose condensates 
studies [5, 6]. This physical situation is also the picture underlying the freeze-out hypothesis used to describe data 
in heavy ion collisions (see section IV A 3). At a given time the system is supposed to occupy a given volume and 
to behave thereafter like an ideal gas of fragments or particles. At variance with experiments in traps, for collision 
experiments the minimum biased distribution (46) at time to is only an ansatz, implying that other observables might 
be needed to correctly describe the density matrix at the initial time to- This more complicated case is studied in the 
next section. 

What is common to both these situations, is that the state (46) will evolve in time. So let us apply the formalism 

developed in section III D assuming that < i?^ > is known at a time while the system is observed at a later time 
t. Since energy is a trivial constant of any Hamiltonian motion, {H, H} = 0, we do not have to introduce a specific 
time for the observation of the energy. 

According to section (HID) to take into account the time evolution, we must introduce additional constraining 
observables 

B^j^^ = {H, = - 1^ ^ {Pn ■ fn + fn ' Pn) 

n 

n 

Since {H,B^^^} — 0, all the other b'"^^ with p > 2 are zero. The above relations are valid both in quantum and 
classical mechanics. Then the most general density matrix corresponding to the time-dependent Gibbs ensemble is 
given by 

D0M it) = ^ E ^ - ^O^r. + ^(pn- ?n + K ' Pn) , (47) 

with 

/?' (t) = /3 + 2Ao {t - tof /m (48) 

and 



vo [t) = 2X0 {t - to) /to. (49) 

Let us explicitly show that the density matrix (47) can be interpreted as a radially expanding ideal gas. Indeed the 
distribution can be written as 

1 ^ , (pn-mho{t)rn) , ^ 
D0,xo it) = ^expJ2 it) ^ - it) ?n (50) 

where the Hubblian factor reads 

_ VQ jt) _ 2Ao {t - to) 
°~ P'{t)~ 0m + 2X{t-to)' ^ ' 
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while the confining Lagrange multipHer is transformed into 



(t) Ao/?m 



A' (t) = Ao - = — 2 (52) 

In the density matrix (50) the term mho (t) correcting the momentum can be interpreted as a collective motion 
produced by a radial velocity /iq (i)r„. This proportionality of the velocity with f„ shows that the motion is akin 
to a self-similar Hubble expansion. As a consequence, when this collective motion is subtracted from the particle 
momentum, the density matrix (50) corresponds at any time to a standard equilibrium (46) in the local rest frame. 

Computing the time evolution of the distribution (50) under the action of the Hamiltonian (45), it is easy to verify 
that the information theory ansatz (50) is, at every time t, the exact solution of an ideal gas initially trapped in an 
harmonic oscillator potential up to the time t{), and then frcx^ly expanding in the vacuum after the confining potential 
has been suppressed. In this case the infinite information wliich is a priori needed to follow the large-time evolution 
of the density matrix according to eq.(21), reduces to the three observables P , p^, f-p + p ■ f. Indeed these operators 
form a closed Lie algebra containing the Hamiltonian operator, implying that the exact evolution of (50) preserves it 
algebraic structure. 

An easy way to follow this time dependence is to compute the evolution of the 3 averages <r^ > , < Pn-fn+fn -Pn > 
and < p^> . Using dt < A >= —{H, A} we can easily show that 

dt <r^>= — <Pn-rn+rn-Pn> 
m 

dt < Pn ■ Tn + rn ■ Pn >= — < Pn > ■ 

m 

The energy being constant, < ^ > is also constant so that the radial momentum < Pn-Tn + ^n-Pn > linearly depends 
upon time while the motion of < f ^ > is uniformly accelerated. This evolution of the constraints corresponds the 
dynamics of /?', and A' deduced above. 

The description of the time evolution when considering unconfined finite systems has introduced a new phenomenon: 
the expansion. One should then consider a more general equilibrium of a finite-size expanding finite-systems with /?', 
ho and Ag as free parameters. Then, if the observed minimum biased distribution at time t is coming from a confined 
system at time to, the three parameters (3' , ho and Ag should be linked to the time to, the initial temperature f3~^ 
and the initial Ag by equations (48), (51) and (52) [25]. 

The important consequence of that is that radial flow is a necessary ingredient of any statistical description of 
unconfined finite systems in the presence of a continuum: the static (canonical or microcanonical) Gibbs ansatz in a 
confining box which is often employed[23] misses this crucial point. On the other hand, if a radial flow is observed in 
the experimental data, the formalism we have developed allows to associate this flow observation to a distribution at 
a former time when flow was absent. This initial distribution corresponds to a standard static Gibbs equilibrium in a 
confining harmonic potential, i.e. to an isobar ensemble. We will see in the following that more complex non-Hubblian 
flows can be associated with additional constraints on the matter distribution. 



C. Deformed unconfined systems 



The systems formed in many experimental situations, ranging from asymmetric atomic traps [6] to non-central or 
incompletely damped nuclear collisions [26], present an explicit deformation. The simplest one is of quadrupolar nature, 
and corresponds to a non-zero < Q2Q > where Q20 is the standard quadrupole moment Q20 = X^n ^zOn where the 
single particle operator reads g20„ = 2z^ — (5:^ + y^) . Therefore, the constraining observable Q20 associated with the 
Lagrange multiplier A2 should be introduced in the definition of the out of (spherical) equilibrium statistical ensemble. 
In addition, following the procedure explained above to take into account the influence of the time dependence, this 
observable should be complemented by its multiple commutators with H 

^Q^ = {H,Q20} = 2pz„Zn - Px„Xn - Pyjn + {p ^ v) (53) 

n 

which is nothing but an asymmetric flow, and 

= {H, S« } = ^ ^ 2pl - pI - pI , (54) 

n 
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a deformation in p space. Again, since {H, Bq^} = 0, all the other Bq^ with p> 2 are zero. 

Then the minimum biased density matrix characterized at a time to by the compactness < > associated with a 
Lagrange multipher Aq, a quadrupole deformation < Q20 > at a time <2 imposed by a second Lagrange multipher A2, 

and by an average energy < P^/2m > , is given by 

D0M it) = exp -P'z it) ^ - A.^^ + ^ iPzJn+ZnPz„) ■ 

exp Y: it) ^ - A^rl„ + ^ • f'x„ + f'^„ • P^^„) , (55) 

n 

where we have introduced the momentum (p±„) and coordinate (f±^) perpendicular to the deformation axis z. The 
distribution at a time t (55) is then function of the spatial deformation parameters = A0 + 2A2 and Aj_ = Aq — A2. We 
can see that the kinetic energy is not determined by a unique temperature but by two different Lagrange parameters 

/3^(t) = P + 2{Xo{t-tof + 2X2(1 -t2f)/m (56) 
/?l(t) = P + 2{Xoit-tof -X2it-t2f)/m (57) 

Finally the deformed expansion parameters read 

iy,{t) = 2{Xoit-to) + 2X2{t~t2))/m (58) 
iy±{t) = 2{Xo{t-to)- X2{t-t2))/m. (59) 

The minimum biased density matrix (55) can also be written as a generalized equilibrium in the local rest frame 

n 1 fl/ / ,N (P2„ - TO^z (i) ^n)^ w .2 /(;n\ 

Di3,XoMit) = y exp2^-/3^(f) — X^z^ (60) 

exp Y it) ^ ^ ^ - Alfl„ (61) 



where the Hubblian flow factors read 



P'At)' P'Aty ^ ^ 

while the confining Lagrange multipliers are transformed into 

^^^'^-^'-^^y ^^^'^-^^-^WJfy ^^^^ 

Computing the exact time evolution of the density matrix (60), and using again the fact that for all k, representing 
both the particle and the axis labels, p^, ff. and [p^fk + VkPk) form a closed Lie algebra, it is easy to demonstrate 
that the density matrix (60) is, at every time i, the exact solution of the dynamical evolution. For example, an 
ideal gas initially trapped in a quadrupoly deformed harmonic oscillator potential up to the time io = ^2, and then 
freely expanding in the vacuum up to the time t exactly follows (60). The time t2 start to play a role if the trapping 
potentials in different directions are not switched-off at the same time. For such a case the density (60) remains the 
exact solution. 

It is interesting to note that the A's are inversely proportional to the size of the system in coordinate space, and so 
are the produced flows. Thus the deformation in p space is larger in the narrower direction in the coordinate space. 

Taking advantage of the commutator algebra of various functionals of the operators p and r this discussion can be 
generalized to more complex multipolar deformations and flows (see appendix A). 



D. Role of an external potential 



In many physical situations the Hamiltonian of the system does not reduce to a simple ideal gas, but also may contain 
explicitly a potential term (j = U ^f„^ , either coming from the interaction with some external fleld, or from a 
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self interaction treated at the mean-field level. In such a case the Hamiltonian becomes H = JZn.Pn/'^''^ + ^ (^n) • 
Considering the additional constraint B!^ it is easy to show that B^^^ = — J2n m ' + " Pn^ remains unchanged 
since |c/ ^f„^ ,f^| =0. The difference starts at the second order, since in addition to the term X^„2^/m^ , 



contains 

2n 



m 



For short time evolutions (up to the second order in 6t = t — to) the minimum biased density matrix corresponds to 
a self similar flow in a modified potential 



1 ^ \ {Pn- niho {t) fn) 
D0M (i) = ^ exp -/?' [t) ^ 5- ^ + U' 



where the time dependent temperature and the Hubblian factor ho are given by eqs.(48),(51), and the effective 
potential tj' is given by 

In the special case of an harmonic potential U ^f^ = the solution can be exactly worked out at any arbitrary 

time. Indeed in this case the B(p'> operators read for p > 1 



B{2p) ^ ^(-l)P(2a;) 

n 

^(2p+i) ^ _^(_i)P(2a;) 



2 2mk 



2p Pn ' ~r fn ' Pn 



m 

n 

where uJ^ = k/m. The distribution at time t reads 

Dp',KAt) = exp -P'it) it) + ^ {pn ■ fn + fn ' k) , 

where the temperature is oscillating in time 

a) = p- Aq (cos 2a; jt-to) - 1) 

as well as the effective constraining field 

,2 



and the coUective radial velocity 

sm2uj{t-to) 
fit) = Ao -. 

Similar to the case of an ideal gas section IV B, this density matrix can be interpreted as a Gibbs equilibrium in the 
rest frame of a breathing system 

1 ^ (Pn- mh' {t) fn) . 

Dp',y (i) = ^ exp ^ -/?' {t) A L- - A' {t) (65) 



where the Hubblian factor reads 



and the effective pressure is given by 
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^ ' I3'{t) rriLO ^ ' 



\'m Ag (sin2a;(t-to))^ 



Due to the confining harmonic potential the distribution is periodic in time. The collective flow is oscillating but 
it is interesting to note that the velocity is at all times proportional to the radius. The expansion deviates from 
self similarity only in presence of an anharmonic potential or of finite range two body interactions, as it is shown in 
appendix B. 



V. CONCLUSION 



In this paper wc have introduced an extension of Gibbs ensembles to accoimt for time dependent constraints, 
in classical as well as in quantum mechanics. The formalism is developed for generic equations of motion and is 
illustrated on both Hamiltonian and dissipative time evolutions of the density matrix. We show that the time 
dependence imposes the introduction of new constraining observables. In the case of an Hamiltonian dynamics theses 
new relevant informations are the multiple commutators of the initial observable with H. This leads to time odd 
constraints wliic;li can be interpreted as collective flows. 

This formalism gives a statistical description of a system characterized by some relevant observables deflned at a 
time at which the entropy has not reached its saturating value yet, as it may be the case in intermediate energy 
heavy ion reactions [26]. Another physical application concerns systems for which the relevant observables pertain to 
different times, as in high energy nuclear collisions where the kinetic energy seems to be still dissipated when the 
chemistry of the system is fixed[9] (i.e. systems presenting difi^erent freeze-out times for diflcerent observables). 

An important result is that any statistical description of an unbound finite system must necessarily contain a local 
collective velocity term. Indeed the knowledge of the average spatial extension of the system at a given time, naturally 
produces a flow constraint at any successive times. This describes the dynamical expansion of the system. This is 
important for the transient unconfined systems formed in coUisional processes which freely expand in vacuum. It 
may also describe out of equilibrium systems bound in a self consistent mean field or confined in traps. Conversely 
a collective flow measurement at a given time can be translated into an information on the system extension at a 
former time. 

In the general case the time dependence of the density matrix can be accounted by introducing an infinite set of 

time even as well as time odd constraints and the relevant information (i.e. the number of independent relevant 
observables or state variables) rapidly increases with time. In some cases however the state variables form a closed Lie 
algebra, and the knowledge of a finite number of observables at a given time allows to predict the total density matrix 
at any successive time. This is notably the case of a system of particles interacting through polynomial one-body 
or two-body potentials. The simplest case is a system of particles or clusters thermalizcd in an harmonic oscillator 
which is suddenly taken off, and freely evolving in the vacuum afterwards. In this case the system experiences a 
self similar expansion at any successive time, and the exact time dependent density matrix is given by an isobar 
canonical Gibbs equilibrium in the local rest frame. This schematic example may have some relevance in the study of 
low energy nuclear collisions, where the so called freeze out hypothesis [23] describes the strongly interacting diluted 
nuclear system as an ideal gas of clusters [2 7]. The presented formalism provide a systematique treatment of the 
confinement and of the successive fiow. 
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VI. APPENDIX A: HIGHER MULTIPOLES DEFORMATIONS 



The reasoning of section IV G can be extended to higher multipoles of rank k, typically considering the operators 
B-"' = J2n'''in^ where fi,„ is the i*'' component of the coordinate of the particle n. Then to have a complete 
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information on the system at any time we have to introduce k additional observables Bf\l = l,...,k given by 

(-1)' kl 



Ml) = y-^) V V 



(2m)' ik-iy.^^^slil-s)l 

This produces non Hubblian flows and non Maxwellian local momentum distributions. For example a B^^^ = 
tensorial deformation, which is part of an hexadecapole constraint, introduces the time dependent constraints sjc''^ = 

-'^Y.n^Px^x'i + xlp^J/m , = ■iT.niPlr.^l + '^Pxr.xlPxr, + xIpIJ/w? , B^x^ = -3E„(Px„in + ^pI„X„Px^ + 

•^Px^XnPx^ + XnP^^)lm^ , B^^ = 24^^jP^/m^. The first term can be interpreted as a cubic flow meaning that the 
expansion is non self-similar. The second term can be seen as a local reorganization of the temperature while the last 
one is an explicit deviation from the Gaussian distribution of velocities i.e. from a Maxwell distribution as expected 
from the Boltzmann factor of an ideal gas. The generalization to several dimensions, i.e. to constraints of the type 
^(0) _ y^" is tedious but straightforward. The B^'^ are linear combinations of products of , p'y and 

with a;* , and z'^ , with i, j, k = 0, . . . ,1 and k^ = i' + i, kj = j' + j and = k' + k. 

It is interesting to notice that considering a generic constraint = Xln/(^")' ^^^^ correction in time is 
B^^^ = X)„Pn ■ ^fi'f^n) + V/(f„) -pn SO that for short time fluctuations {t — to) = 5t the statistical ensemble 
reads (at the first order in 5t) 

2 



with 



1 I (Pn~^(fn)) 

D,At) = ^ exp 5: ^ + U 




A 



(4) = StVU{ir^) 



which is akin to an equilibrium of particles in the external scalar and vector field U and A. 



VII. APPENDIX B: REAL GAS 

Let us now study the case of a real gas with H = K + V with V = X^nn' ^ two body interaction depending 



only upon the relative distance f„ 



Let us investigate the minimum biased distribution including 



a compactness observable known at a time tf/. V depending only upon f , sj^-* remains unchanged: B)^ 
{K, R^} = — ^ X^n^" ' + • Pn while B^j^^ contains an additional term 

^ m m 

nn nn' 

and thus reads 



n nn' 

In case of an harmonic interaction the operators only contain quadratic terms J2nn' ^nn' ^^'^ Snn' ' 

Pnn' , with pnn' = Pn — Pn' ■ lu this casc the time evolution can be taken into account by a suitable time dependent 
reorganization of the temperature and the introduction of a time odd constraint, the radial flow. However, for any 
other interaction B^^ modifies not only the temperature but also the two-body interaction. If we define V'{f) = 
V{f) + rdrV{f)/4: and we work out the third order term: 

2 



SL') = 4. E {nfnn'),pi.} - ^ E {pi^fnn'drV{fnn')} 



nn 



2 / . Vnn' ' "^nn' ^ ~r ^ 'fnn' ' Vnn' 

^ , '^nn' "^nn' 

nn' 
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we can sec that the time dependence of the process induces an effective momentum dependent two-body interaction. 
An interesting phenomenon occurs at the level of the next order. Indeed while the kinetic energy term in the 
Hamiltonian leads to a first term in = {H, -B^^} 



^i"' I I Inn' I nn' i 

nn' \ ) 

which is again a akin to a two body interaction, the interaction part produces a three-body term 
{^,^3)^ = \ ^ \vir^,VPn-^r.n'^^^^^'^^ 

nn n v J 

= ^ X] drV{r^^'')drV'{rnn') 
nn'n" 

n.r) ' n." ^ 



VIII. APPENDIX C: DETAILS ABOUT THE RPA 



In eq. (30) we have expressed the RPA matrix as H = W — AfV where the self consistent mean field is defined as 
W II • >=|| [W, ■] >, the residual interaction V || • >=|| > and the density operator ^f \\ ■ >= jj [p, •] >. 

To show the hermiticity of these operators one should go back to the definition <C A || B+ || C ^= 
(< C II ^ II ^ >)* • In the case of W and , using {TrC+ [B, A])* = TrA [B, C] valid when B is hcrmitian B+ = B, 
we get <SC >1 II S+ II C 3>=<C A || S || C . Concerning the residual interaction, using the hermiticity of W and its 
definition we can easily get an explicitly hermitian form V[i^2] = (^'^ ^ / 9 Py^d pj^y 

Then we have expressed the generalized constraints of the time dependent RPA problem in the small amplitude 
limit as 



Ai = a^fi^,+alfi+^ (68) 
B(, = iw^^a^^C^^ - iuj^^al^C+ (69) 

and we have interpreted A^ as a collective coordinate associated with the creation of a collective mode through the 
(7+ operator, and Bi as the associated momentum. 
Let us justify this interpretation. 

If is an eigenvalue of the RPA matrix corrcponding to the cigcnstatc 6pi,, it is easy to show that 5/5+ is solution 
of CO* II (5p+ ;»= —H II ^ so that || ^ is also an eigenstate of the RPA matrix associated to —ui*. Thus RPA 
solutions can be grouped by pairs. 

Then, because of the commutator structure of the dynamical equation, the RPA equation do not propagates the 
diagonal terms dpu, so we can focus on the off diagonal terms . Thus we can introduce a collective operator Q such 
that II 5p »=|| [p, Q] A/^ II Q » . In fact the density variation dp can be interpreted as produced by a unitary 
transformation generated by the operators Q and Q~^: p — »■ e~*'^('^+Q"'")pg»-\(Q+Q''") = p + iX[[p,Q] + [p,Q~^]). Thus 
the relation holds lo^, \\ >= N'^TiN || Q,. > . 

If \W, p] = 0, then Af~^HN' = W — VAf ='H~^ meaning that || Qi/ » is an eigen vector of W+ associated with the 
eigenvalue w^.This means that, when [W, p] = 0, cj* is also an eigenvalue of the RPA matrix and, since the RPA 
eigenvalues appears by pairs —uji, is also a solution. Requiring the equality of these to solutions with the pair Ui, and 
-to* shows that u)i, is real or purely imaginary and is associated with a second solution — w^. 

If we now introduce the dual basis <C C„ || such that <C C„/ || dpi, ^= J^/i/ and the associated closure relation 
J2„ II Sp^ X C„ 11= 1, the relation < C,; || 7i < C^, || is easily demonstrated. Then w* || >= 7^+ || C„ > 
also holds showing that, when [W, p] =0, = Qv if is real. 

This shows that the dual of the eigenvector || Spi, S> is the operator responsible of the excitation of the collective 
mode LVi, as indicated in equation (69) with a^,^ =<C A^ \\ Sp^, 
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